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ABSTRACT 

We consider the one-dimensional partially asymmetric exclusion model with open bound- 
aries. The model describes a system of hard-core particles that hop stochastically in both 
directions with different rates. At both boundaries particles are injected and extracted. By 
means of the method of Derrida, Evans, Hakim and Pasquier the stationary probability mea- 
sure can be expressed as a matrix-product state involving two matrices forming a Fock-like 
representation of a general quadratic algebra. We obtain the representations of this algebra, 
which were unknown in the mathematical literature and use the two-dimensional one to derive 
exact expressions for the density profile and correlation functions. Using the correspondence 
between the stochastic model and a quantum spin chain, we obtain exact correlation func- 
tions for a spin-i Heisenberg XXZ chain with non-diagonal boundary terms. Generalizations 
to other reaction-diffusion models are discussed. 
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I. INTRODUCTION 



One dimensional reaction-diffusion processes recently have attracted much attention for a 
variety of reasons. Pure diffusion models have been studied in relation with interface growth 

traffic-flow 0, the dynamics of shocks @,|], and magnetophoresis of tagged polymers 0. 
More general reaction-diffusion models are of interest from a mathematical point of view due 
to their relation to integrable quantum chain hamiltonians 0. It is interesting to note the 
important role played by the boundary conditions in these models 0, which completely control 
the physics in some cases. For the case of two-state models for example, in the corresponding 
quantum chain hamiltonians (which are XXZ models) the boundary conditions generally break 
the particle number U{1) symmetry and are not easily treatable by the usual methods like the 
Bethe Ansatz MM. The problem is that although the chains can be shown to be integrable 



10| , the Bethe Ansatz so far has not been constructed due to the lack of a reference state. 
An important step forward in these types of problems was made by Derrida, Domany and 



Mukamel []TT| in the case of completely asymmetric diffusion with particle injection at one 
end of the chain, and particle extraction at the other end of the chain. They showed that 
there exists a recursion, which relates the probability distribution of the steady state for a 
L sites to the one for L — 1 sites. An equivalent formulation of this property was given by 
Derrida, Evans, Hakim and Pasquier (DEHP) [T^, who demonstrated that the probability 
distribution can be written in a factorized form with coefficients that are not c-numbers but 
(infinite-dimensional) matrices. For the two-state model there are two matrices which form 
a Fock representation of the quadratic algebra. Using representations of this algebra one can 
compute in principle all correlation functions. In particular the density profile was determined 
in [^,0, and in a special case (when the injection rate is equal to the extraction rate) even the 
two-point correlation function was obtained. If one considers the more general problem 
with particle injection and extraction at both ends, and partially asymmetric diffusion, the 
DEHP approach is still applicable, but the representations of the quadratic algebra were 
not known. In a remarkable paper []T5[, Sandow was able to compute some important matrix 
elements in the envelopping algebra, which allowed him to compute the currents and to obtain 
the phase diagram, which coincides with the mean-field predictions. 



In the present paper we start by determining all Fock representations of the most general 
quadratic algebra, which depends on seven parameters (section II). They might be of interest 
in other physical contexts as well. It turns out that the representations can be either finite- 
dimensional or infinite-dimensional. For each finite-dimensional representation one obtains a 
constraint equation for the seven parameters. This constraint depends on the dimension of the 
representation. The matrix elements of the two matrices appearing in the quadratic algebra 
are given by recursion relations. We show that for some special cases these recursions can be 
easily solved. 

Next we review the connection between the steady-state probability distribution and the 
ground state of certain (in general non-hermitian) quantum chains (section HI). In section 
IV we give a summary of the DEHP Ansatz and establish the connection with the quadratic 
algebra discussed in section II. In section V we consider the most general master equation 
for one-dimensional systems with two-body interactions (one has twelve independent rates) 
and particle injection and extraction at both boundaries and apply the DEHP Ansatz. One 
gets a quadratic algebra and two additional quadratic relations on the matrices. This implies 
that only finite-dimensional representations have a chance to be useful. We found that there 
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exists a one-dimensional representation with three conditions for the twelve rates and the four 
parameters describing the injection and extraction of particles at the ends of the chain. The 
question of the existence of higher-dimensional representations and their physical relevance is 
left open. In Appendix B we study the applicability of the DEHP formalism to the steady 
state of a master equation with three-body interactions. We show that is this case one obtains, 
as expected, cubic algebras. Their representations and physical interest remains to be studied. 

After these mathematical investigations we turn to a detailed study of the problem of 
partially asymmetric diffusion with particle injection/extraction at both boundaries. We 
start with a review of the known results in section VI, using the phase diagram obtained 
by Sandow |]l5l as a basis. In the following sections we concentrate on the application of the 
two-dimensional representation of the quadratic algebra to concrete calculations. In sections 
VII and VIII we show that as a result of the constraint equation for the existence of the rep- 
resentation one can cover parts of the phases Aji, Bjj, the complete phases Aj and Bj as well 
as the coexistence line (in terms of the definitions of Schiitz and Domany [l^). The calcu- 
lation of the density profiles and two-point correlation functions in the low- and high-density 
phases is presented in section IX. One remarkable result is that the density around the center 
of the chain has a simple expression in terms of the parameters of the problem and that it 
coincides with the mean-field results (which are derived in Appendix F). The density profile 
and two-point function on the coexistence line are presented in section X. 

Some by-products of our investigations are presented in the appendices. In Appendix C 
we give some identities concerning normal-ordered expressions of q-oscillators. In Appendix 
D we show how the DEHP Ansatz can be used to construct irreducible representations of the 
quantum group Uq{SU{2)). We close with a discussion of our results and some remarks on 
the time-dependence of correlation functions. 



II. FOCK REPRESENTATIONS OF THE QUADRATIC ALGEBRA 

We are interested in Fock representations of the most general quadratic algebra 

xiA^ + X2AB + X3BA + x^B'^ = X5A + xqB + X7 , (1) 
A\V)=0, {W\B = 0, {W\V)^0. (2) 

Here Xi are complex parameters and quantities of physical interest are given by vacuum average 
values of monomials written in terms of A and B, e.g. 

{W\A'~'B'\..B^-\V) . (3) 

Obviously (|^) generalizes the algebra of creation and annihilation operators and its q- 
deformations. As far as we know the problem formulated above has not been considered 
in the mathematical literature (presumably because up to now there was no motivation to 
do so). General quadratic algebras were studied in but no Fock- representations were 



considered. We will show in the present work that solving the above problem allows for the 
computation of concentration profiles and various correlation functions in the physical problem 
of partially asymmetric diffusion with open boundaries. 

One can ask the question about the conditions on the Xj's in (|^) such that algebra de- 
termines the vacuum expectation value (^. This implies that the system of equations for 
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words of length two, three, etc have solutions. Direct calculations show that an infinite set of 
inequalities have to be satisfied 



X2 7^ , X1X4 — X2 7^ , X2{xl — X1X4) + (x3 — X2)xiX4 7^ , etc . (4) 



In section IV we will give, in a different parametrization, a simpler expression for these condi- 
tions. Instead of solving linear equations for words of different lengths it is useful to look for 
matrix representations of the algebra. Once those are known the calculation of the quantities 
d^) is simple. We are interested for obvious reasons in the representations of the smallest 
dimension because this is sufficient to compute the relevant quantities. As we are going to 
show the representations of the quadratic algebra are infinite-dimensional unless there exist 
certain constraints on the parameters Xj. The simplest such constraint is X7 = 0, for which 
the representation is one-dimensional: A = B = 0. From now on we will take x-j 7^ 0. First 
we consider the case X5 7^ 7^ xg. We then define 



A=^A 
X7 



B = ^B 
X7 



in terms of which the algebra reads 

ziA^ + Z2AB + z^BA + z^B"^ = A + B + l 
where the Zi are given in terms of the Xj and where 

{W\B = A\V) = 0, {W\V) ^ . 



(5) 



(6) 



(7) 



It is convenient to define 

XiXe 



Z2 



X2X5 



Z2 X2 Z2 X2XQ 



Z2 



X^Xq 



(8) 



One can show (see Appendix A) that through a similarity transformation the matrices A and 

/1\ 



B can be brought to a tri diagonal form with |V^) 



this fact together with equation we obtain 
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The quantities a„, a^, hn and /„ are given recursively. First the a„'s are to be determined from 

1 + AOn-l 



Or 



1 - rj^Un-l 



, ai = 



(11) 



Next one determines a„ and 6„ from 



M 



A + r]^an r]{l - an) 
1 

ai = = 61 . 



M 



a„ \ W 1 + 

, &n / ^2 \ 1 + r/ar 
-77(1 + Aa„) -l+?7^a„ 

-1 + ?7^a„ -^(1 + Aa„) 



(12) 



Finally the /„'s are then given as 

J n J n — 1 



1 - 2?7^a„ - AT^^a^ 
/o = 0. 

For later use we give the first few values 

1 + A 



+ 



1 - 2?7^a„ - Xriial 



(13) 



"2 = 1, "3 
02 



1-r/e ' 
1 + r/ , l+i 



^2(1 - ^0 ' ^2(1 - ^0 



fl — 5/2 
Z2 



A + 1 



zi(l-2r/^-Ar70 



2:2 + 



:i+o(i+^) 
(1 - riif 



(14) 



We note that for A = —1 it follows that /2 = and we thus obtain only two-dimensional 
representations like for one fermion (observe the appearance of an anticommutator in (^). 
The recursion (|rT|) for a„ can be solved by redefining 



(15) 



where the m„'s satisfy the following recursion relation 
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1 + A 

Un+l + Un-1 — / .. , ^ Un , 



Ml 



1 , U2 



(16) 



This is recognized as the special case X 
nomials Un+i(x) + f/„_i(x) 



l+A 



2i^x+riO recursion relation for Chebyshev poly- 

2xUn{x). Using the representation Un{x) 



and taking into account the initial conditions we arrive at the result 



sin{(n— 1) arccos(x)H 
sin(0) 



1 



sin {{n — 1)6 



arctan 



- (1 - A) 



A 



sin ((n- 2)0 + 

, e = 



arccos 



1 + A 



2VA+^ 



(17) 



The recursion ([T^ ) for a„ and bn can be decoupled into recursion relations for a„ = y/^an + y/fjbn 
and bn = v^ctn - ^/nbn■ 



— 9n+l^n + ) — 9n+l^n + ^n+1 ' 

where /i^+i = ci ± C2, (7^+1 = C3 ± C4, and 

:r]{l + Aa„)(l + ^a„) + (1 - ?7^an)(l + ?7ar 



(18) 



ci = -Ve- 
C2 = -y/v- 



C3 



C4 



;22(r7e(l + Aan)2-(l-r7e«„)2) 

1 + ^Ctn)(l - ^^ftn) + ^(1 + ^an)(l + Act^ 

^2(r/e(l + A«„)2_(i_^^^^)2) 
(A + r;0(l-^e«n) 



r?e(l + Aa„)2 - (1 - r]^any ' 
v/^(A + l)(l + (A-l)a„ + r7ea2) 



r/e(l + Aa„)2-(l-r^ea„)2 



(19) 



It is hard to simplify the recursion relations further. Using the expression ([17D for the a„'s and 
([T9| ) one can derive formulae for a^, 6„ and hence /„. The resulting expressions are obviously 
very cumbersome. From the expressions (^ and (|I0|) for A and B it follows that the condition 
for having an n-dimensional representation is simply = 0. As one can see from the form of 
the recurrence relations this constraint is a complicated function of A, ^, rj and Z2- This is the 
reason why we will use for applications only the two-dimensional representation (/2 = 0), for 
which the matrix elements of A and B are given by ([T^). 



Let us now consider the case X5 = Xg = 0. The cases where only X5 or xq vanishes can be 
studied in a similar way and thus will not be considered in detail here. The algebra for case 
X5 = Xq = is 

ziA^ + Z2AB + Z3BA + z^B^ = 1 , (20) 
where Zi = f-. We define, like in (||) 



5 



Z2 



A = - 



£3 

Z2 



v = — 



£4 

Z2 



(21) 



The infinite-dimensional representation of (pOf ) is of the form (^-(|13) with vanishing diagonal 
terms a„ = 6„ = , Vn, where r],^, X are defined in (^TD and where q;„ and fn are given by 



Or 



1 + \an~i 
1 - r]^an~i 

fn-1 



«! = 



2 A + 2r7^a„„i 



+ 



1 



1 - 2?7^a„ - Ar^^a^ ^2(1 - 2r]^an - Xr]^al) 







(22) 



Apart from this infinite-dimensional representation there are two kinds of finite-dimensional 
ones. The first kind is simply obtained by imposing the constraint /at = on the parameters 
Zi, which leads to the decoupling of an x block in the upper left corner of the infinite- 
dimensional representation of A and B discussed above. The resulting finite-dimensional 
representation is given in terms of N x N matrices A and B with vanishing diagonal elements. 
The matrices A and B of the second type of N x N representation take the following form 
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where a„ and /„ are determined by the recursion (^). The representation (]23|), ( p4D exists 
provided that 

c^nV^ = 1 • (25) 
The variables a^r and bj^i are obtained from the equations 

(-A - 2^?7aAr_i + r]^a%_^) /^.^ = + r]b% + (A - l)aNbN + — ■ (26) 

We are going to close this section with two cases, for which the recurrence relations can 
be solved in a trivial way. 

• If A = — 7^—1 (this case is as we will see physically interesting) the following simpli- 
fications take place for X5 7^ 7^ Xe: 

1 + 27] - A^ 1+r/ 
an = a = — — —-- , Vn > 3 , ai = 0, 02 



z,{i + xy - - ^2(1 + A) 

bn = b = ^^^^^^^^ , Vn>3,6i=0, 62 ^ ^ 



^2(1 + A)2 ' - ' ^ ' ^ ;22(1 + A) 

2 r2 



In- J - n I , vn ^ d 



(1 + A)^ 



Note that in this case there exist only 2-d (/2 = 0), 3-d (/s = 0), and infinite-dimensional 
representations ! If X5 = Xg = 0, one has an infinite-dimensional representation with 
vanishing diagonal elements a„ = 6„ = , Vn with 

«„ = lVn>2. /? = i./J = -^./„^ = -^.V„>3. (28) 
where ?7, ^, A are given by (^TJ). 

o = ^ = e 

Here it is possible to choose B = A\ and 

where we have assumed Z5 7^ 7^ xg. A similar simplification holds if X5 = or xq = 0. 
Of particular physical interest is the case A = 1, which exhibits additional simplifications 

n — 1 2 ^ ^(^ ^1) 
= n - 1 , a„ = , /„ = \ 2 • (30) 

Z2 Z2 Z2 
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III. THE MASTER EQUATION AND THE QUANTUM CHAIN HAMILTONIAN 



Let us consider a one-dimensional open chain with L sites. On each site k{k = 1,2, ...L) 
we allow for two configurations described by means of the variable r^, which takes the two 
values and 1. For Tk — the site k is empty (vacancy), for Tk — 1 the site k is occupied by 

a molecule A. 

At time t the probability to find a certain configuration of molecules and vacancies on the 
chain is given by the probability distribution 

PL{Tl,T2,...TL\t) . (31) 

If we assume that interaction between molecules is described by two-body processes only 
(thrcc-body processes are considered in Appendix B), the time evolution of the system is 
given by a master equation of the form 

Qp L-l 

-Kf- = - I] (^fc,fe+l)rlrfc+i ^Liru ^2 ■ ■ ■ Tfe-l, 7ik, 7ik+l, ■ ■ ■ Tiji) 

^ fc=17fc,7fc+i 

- E (^i)n Pl{1i, r2...r,\t)-Y. {hL)T, Pdn, . . . t,.„ n\t) , (32) 

71 IL 

where the boundary contributions hi and Hl describe injection (extraction) of particles with 
rates a and 5 (7 and /3) at sites 1 and L 





] ^ \ _ ] , (33) 

and where 

^ -[^^Ml?rt:i else 

Here [Ffe^fe+iJ^^'^^+j represents the probability per unit time that the configuration (7^,7^+1) 
on neighbouring sites k and k + 1 changes into the configuration (r^, Tk+i) and Y^' denotes the 
sum where the term (7^,7^+1) = {jSk, l3k+i) is excluded. The following processes are included 
in the master equation (with a vacancy (r = 0) and A a molecule (r = 1)) 
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Diffusion to the right: A + - 

Diffusion to the left: + A- 

Coagulation at the right: A + A ■ 

Coagulation at the left: A + A ■ 

Decoagulation at the right: A + - 

Decoagulation at the left: + A - 

Birth at the right: + - 

Birth at the left: + - 

Death at the right: + A ■ 

Death at the left: A + 0- 

Pair-annihilation: A + A 

Pair-creation: + - 



+ A (rateP 
A + (P 
+ A 



^ A + 
'A + A 
'A + A 
+ A 
A + 
^0 + 
* + 
■^0 + 
A + A 



(P" 



01 



(r 



(r 



11 



(poo 



01 



(r 



(roi 



00 



(r 



(p" 



00 



(r 



(35) 



React ion- diffusion models of the type described above can be mapped to quantum spin chains 
in the following way 0: a basis of the quantum-mechanical Hilbert space Ti. (isomorphic to 
the tensor product ®^^|C^) is defined as 



\{r}) = \n...r,) , 



(36) 



and the inner product is taken as 



This induces a map of the probability distribution to a state in Ti 



and the master equation 



{r} 

^) then implies an imaginary-time Schrodinger equation 

d\P) 



(37) 



(38) 



dt 



-H\P) . 



(39) 



Here if is a quantum hamiltonian defined in terms of a basis E'j^^ (which can be represented 
as 2 X 2 matrices with entries E°'^^s = ^a-y^ps) of quantum operators on the fc'th site of the 
lattice via 



L-l 



k=l 



(40) 
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where hi, hi and i?fc,fc+i are defined in (|34D . Note that in general the hamiltonian H will be 
non-hermitian. It is easy to see that 

(0|^ = 0, (41) 

where (0 | is given by 

(o|=E(MI= ((1)® ' ^^2) 

Using (|39|) it follows from ( ^T]) that (0| is a left "stationary" state. Assuming that this is the 
unique left "stationary" state and given a unique right stationary state 

|o) = E^^(M)IM) (43) 

{r} 

the average of the observable X (ri, . . . tl) = X({r}) is defined as 

W = E^(M)^.(W) 

{r} 

= (0 I X I 0) . (44) 

An example, with which we will be concerned in most of this paper is the case of partially 
asymmetric diffusion, which corresponds to the choice of rates 

[r,,fc+i]J? = p, [rk,k+i]fo = (i (45) 

(all other rates are taken to be zero) the quantum hamiltonian H obtained by the above 
mapping is related to an XXZ spin chain by a similarity transformation 

Hxxz = UHU-' , f/ = n + Ej'AQ^-') = U i (46) 

where Q = y^, A is a free parameter and 

1 _ _ 1 

-^Hxxz = -2 E 2(^/^7+1 + ^7^;+!) + 2(2 + 

pq ^ j=l L ^ 



^-^(^1+1 -^D- ^(2 + 2"') +Bi + Bl 



_ Aa + 7 a + 7 

-Di = cTi cTi aT- 1 

2y/pq y/pq Ay/pq 2y/pq 



Bl = (J^——- a^^=Q -a^——Q + — — . (47) 
2Wpq y/pq A^/pq 2^/pq 



This is the Uq{SU{2)) invariant quantum spin chain [|T7| with added boundary terms Bi 
and Bl. Notice that the boundary terms contain nondiagonal contributions {af, a^) with L- 
dependent coefficients. In the absence of the boundary terms the spectrum of the hamiltonian 
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is massive. As is shown below the boundary terms with generate phase transitions with 



massless phases. Although the hamiltonian (|T7|) can be shown to be integrable [10|, the Bethe 
Ansatz so far has not been constructed due to the lack of a reference state. 

We note that the similarity transformation ( |^ ) does not change averages of observables 

(0|X|0) = (0|f/-if/Xf/-if/|0) = u{0\Xu\0)u . (48) 

Thus zero-temperature equal time correlation functions of the XXZ quantum spin chain and 
(stationary-state) averages of the partially asymmetric diffusion model are related in the fol- 
lowing way 

{a]) = u{0\a]\0)u = 2{t,) - 1 , 
{a]al) = u{0\a^al\0)u = 4(r,-r,) + 1 - 2(r,) - 2{n) , 

^^j^klconn ^irj^k) conn i (49) 

where conn denotes connected correlation functions. This means that all results concerning 
averages in the partially asymmetric diffusion model obtained in this paper can be immediately 
applied to the case of the XXZ chain described above. 



IV. THE DEHP ANSATZ 

In a remarkable paper it was shown that for the case of asymmetric diffusion the 
problem of determining the probability distribution Pl{ti, . . . t^) for a stationary state can be 
formulated in a completely algebraic framework. We now briefly review the relevant results. 
All rates except [Ffc^fc+ijQ? = p and [Ffc^fc+ij^Q = q {k = 1, ... L — 1) are taken to be zero, and 
the boundary conditions are chosen according to (|33|): particles are injected at sites 1 and L 
with rates a and 5 and extracted with rates 7 and (3 respectively. The algebraization of the 
problem of determining the unnormalized probability distribution Pl{ti . . . r^) of a stationary 
state is performed in two steps: one first makes an Ansatz Pl{ti, . . .t^) in the form of a 
matrix-product state [p!^ , p!^ 



PLin, tl) = {W\ YlinD + (1 - n)E)\V) . (50) 



i=l 



Here D and E are in general infinite-dimensional matrices and {W\ and \V) are vectors 
connected with the boundary conditions. The normalization factor is obviously given by 

Zl = {WIC'-IV) , C = D + E. (51) 

In the second step the following sufficient conditions for Pl to be a stationary solution of the 
master equation are imposed 

(^fc,fc+l)r',Tr+' ^l{ti, T2... Tk-l, 7fc, 7fc+l, n+2 ■ ■ ■ Tl) 

XrkPL-l{ri . . . Tk-lTk+l ...Tl)- Xr^^^PL-lin . . . TkTk+2 . . . Tl) , (52) 

J2ih)^rlPLillT2 ...Tl) = -Xr,PL-l{T2 • • • Tl) 
71 

Y,{hL)llPL{T^ ■ ■ ■ TL-11l) = Xr.PL^M ■ ■ ■ ^L^l) ■ (53) 



7fe.7fe+i 



7i 
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Inserting (^) into (|52D , (|53D leads to algebraic relations between the matrices D and E and 
leads to conditions for the action of D and E on \V) and {W\ |]12[ (one finds that 
and then sets xi = 1) 

pDE - qED = D + E 
i/3D - 5E)\V) = \V) 
{W\{aE--fD) = {W\ 

{W\V) ^ . (54) 



It is easy to see that for aP = •jS no representations of (|54D exist. This can be seen by 
considering the inner product 

{W\{aE --fD)\V) . (55) 

Evaluating (|55|) once by acting to the left and once by acting to the right using that aE—'-fD = 
— ^{PD — 6E) (which holds because a/3 = •jS) we obtain 

{W\V) = -^{W\V) , (56) 

which has the only solution = 0. This means that from now on we can constrain 

ourselves without loss of generality to the case a/3 7^ Moreover one can show that all 
vacuum expectation values 

{W\D'''E'^\..D''"\V) (57) 

are determined by (Q) if the following inequalities are satisfied 

/a/?-gV7^0, k = 0,l,2,... (58) 

The proof uses the construction of [jr3| and one can show that if (^) is satisfied a representation 
exists even if its dimension is not the one of the smallest representation. In order to obtain 
Pl it is now necessary to find matrices D and E together with vectors {W\ and \V) obeying 

For later use we note that physical quantities like the current J, density profile {tj) and 



two-point function (TjTk) can be evaluated in the following way |12 

{W\C^-^V) 



J 



{W\C^\V) ' 
{W\C^-^DC^-^\V) 

(w\cW) 



^^^^'^ = wmv) ' ^ ^ 



^Actually there is one exception: if p = q and a(3 = 
representation {D and E are numbers) with xq = 0, D = 

12 



jS there exists a trivial one-dimensional 
IE. 



where C = D + E. 

It is possible to determine certain matrix elements of representations of the six-parametric 
algebra ( |54D directly. This was done by Sandow [|T5| who then was able to determine both the 



current J (the computations of which involves only the matrix C) in the infinite volume limit 
L ^ oo and, remarkably, the phase diagram of the system. In order to compute the density 
profile and the two-point function a much more detailed understanding of the representations 
is needed. In order to study representation theory of ( |5^ we introduce two operators A and 
B, which act trivially on {W\ and \V) respectively 

A = (3D-6E-1 , B = aE --fD -1 , 
{W\B = = A\V), (60) 

If A and B are known, one can get D and E since a/3 7^ 7^. A and B are seen to obey the 
quadratic algebra @ discussed in section II with 

Xi = (j9 — g)a7 , X2 = paP — q^yd , X3 = p^5 — qaP , X4 = (p — q)pS , 
2^5 = (a + 7) («/5 - 7<^) - (p-q) ilia + 5) + a(/5 + 7)] , 
X6 = (/3 + 6){a(3 - j6) - {p - g)[/3(a + 6) + 6{(3 + 7)] , 

X7 = {a + (3 + + 6){a(3 - 7^) -{p-q){a + 6){(3 + 7) . (61) 

Note that the seven parameters Xi are not independent, since they depend only on six 
variables p, q, a, P, 7, 6. Conversely the algebra (0) can be brought to the form (K^ if 



(X2 + xsf - 4xiX4 ^ . (62) 

One can easily verify that the inequalities ( p^ ) and (^) correspond to those of ( [55| ) for k = 
0,1,2,3. 

We now observe that in the case g = of completely asymmetric diffusion we have A = —77^ 
and the matrices A and B (which we recall were defined as A = ^A and B = —B) have 
the simple forms given by (pTj). The origin of this simplification can be traced back to the 
representation theory of quantum groups which also simplifies drastically in the crystal basis 
q = . Another case for which A and B have a simple form is p = q (symmetric diffusion) . 
Here we have = tj = 0,A = 1, and the representation is given by (|30|) . Finally, if a = /5 = p—q, 
7 = 5 = 0, we can define 



I^A , at = J^B (63) 
p — q \ P — q 



and the algebra (0) is rewritten as a Q-oscillator algebra [pT| ,p^,p5 



aa^ — Qa'^a = 1 , (64) 

where Q = |. The vectors \V) and {W\ turn into usual Fock-vacua (0| and |0) defined by 

a|0) = (0|at = 0. Some observations about this case can be found in Appendix C. 

As noticed before in section III, in the absence of boundary terms the hamiltonian (^Tf ) is 
Uq{SU{2)) invariant. The ground state of this hamiltonian is L -|- 1 times degenerate (recall 
that L is the length of the lattice) corresponding to a L -|- 1-dimensional representation of the 
algebra. We demonstrate in Appendix D that this representation can also be found through 
the DEHP-Ansatz. 
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V. FURTHER APPLICATIONS OF THE DEHP ANSATZ 



It is an interesting question, to what extent the Ansatz (^) can be used to describe more 
general reaction-diffusion models of the type (^). Inserting (|50|) into the master equation 
(|53|) for a general reaction-diffusion process defined via (0)-(^) (note that we take all rates 



constant throughout the bulk, i.e.. 



[^k,k+l]'^S 



7<5 / ) 



we obtain the algebra 



n 



ED 
DE 



( \ 

XqD — X\E 
—xqD + xiE 



V 







/ fOO-LFOO I pOO 
/ Ol"'"-'- lO"!"-'- 11 



V 



^ 01 

_roo 

10 



pOO 

11 



■^01 
00 



/ 

pOl 
J- 00 


plO 

00 


, pOl.pOl 


plO 
^ 01 


pOl 
10 


pio.rio J 

00^-*- 01^ 


pOl 

11 


plO 

11 



plO 

J- 11 



_Fll 
00 

^ 01 

10 



(65) 



pii_l_pii_l_pii ; 
00 01 -"^ 10 / 



whereas the boundary conditions (^) impose the following conditions on the vectors \V) and 
{W\ 



Y.ih^):^l{W\[^^D + {l-^^)E] = 

71 

J2{hL)rAlLD + {l-^L)E]\V)=Xr,\V) . 



(66) 



7i 



Here hi and hi are given by (|33|) and thus (^6D are independent from (^). From (|6g) it follows 
that xq = —Xi, so that we can choose xq = — 1, xi = 1 by fixing the overall normalization in 

dH). 

It is easily seen that only three equations of (^) are linearly independent. These can be 
cast in the form 



niDE + K2ED 
KeE^ 



D + E , 
K4DE + K5ED 
K^DE + KsED 



(67) 
(68) 
(69) 



where are given in terms of the rates r°f . These equations can be viewed in the following 
way: (^) is the basic requirement as it cannot be eliminated by adjusting the rates, whereas 
(|68| ) and (|69D are additional relations in the algebra which are absent in the simplest case 
where K3 . . . Kg are chosen to be zero (by adjusting the rates). This simplest case corresponds 
to partially asymmetric diffusion and will be studied in detail in what follows. The important 
point is that the set of representations of (|67|)-(|69|) is a subset of all representations of ( |67D 
for arbitrary ki and K2- This means that in all cases physical quantities can be determined 
by using a representation of (|67|) only, and then impose the further relations on the matrix 
elements of D and E entering the computation. The existence of solutions of the complete 
system (|66D, (|67D-(|69D is established for the simple case of one- dimensional representations of 
in Appendix E. The question of existence of a two-dimensional representation is still open. 
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VI. ASYMMETRIC DIFFUSION: KNOWN RESULTS 



Before we turn to the derivation of our results for current, density profile and correlation 
functions of the partially asymmetric exclusion model we give a short review of some important 
previously known exact results. So far exact results have mainly been derived for the case of 
completely asymmetric diffusion with injection of particles at one boundary, and extraction 
at the other. In our notation this corresponds to the choice = q = 'j = 6, p = l. This 
corresponds to the infinite-dimensional representation given by (P^ with = X = rj = ^. The 
phase diagram for this case is of the form given in Fig.l []TT],|T2],|T3[ . Note that in order to make 
the connection to the partially asymmetric case easier we have plotted the phases as functions 
of = —1 + ^ and = — 1 + instead of a and p. There are three main phases: a 

high-density phase A, a low-density phase B, and a maximal-current phase C. Phases A and 
B are further subdivided into Ai, An and Bi, Bu respectivel y |T3| (see below; note that we 
have changed notations by switching A and B as compared to [13| in order to comply with the 
notation of ||15|). Phases A and B are separated by a line which is called "coexistence line". 




FIG. 1. 

Phase diagram of the completely asymmetric exclusion model. 

The currents in the three phases are given by (L ^ 1) 

Phase A : J = - f3) , 
Phase B : J = a{l — a) , 

Phase C : J= \ . (70) 
4 

15 



The density profile in the center of the chain (j ~ -I, L ^ 1) is of the form 



Phase A (high density) : (tj) 

Phase B (low density) : (tj) 

coexistence line : {tj) 

Phase C : {tj) 



1 + «:+(/?) 
1 



1 + /t+(a 
a + (l 



2a) 



(71) 



The subdivision of phase A into Aj and Ajj (and similarly B in to Bj and Bjj) was proposed 
in 



chain, which for L ^ j 1 is of the form |]12 



and is based on an analysis of the behaviour of the density profile near the ends of the 



Phase Ai : (tj) 
Phase An : (r, ) 

Phase Bi : {'T'L-j) 



1 + /€+(/?) 



Phase Bii 



Phase C 




2^ J 



(72) 



Note that if /3 = | in phase C there is no j-dependent correction terms in the density. The 
mixed notation in terms of K+(a), k+(/3) and a, (3 has been chosen deliberately and is based 
on universal behaviour in the partially asymmetric case (see below). 

Correlation functions for the completely asymmetric case (0 = g = 7 = 5, p = 1) and 
a = 1 = (3 were obtained in [Q. This corresponds to the point = n+{a) = k+(/3) in the 
phase diagram (this is the Q-oscillator representation ( |63D with Q = 0). In the thermodynamic 
limit L ^ 00, ki 1, k2 ^ I, kj fixed, the connected two-point function was found to exhibit 
an algebraic decay 



1 



4(A;iA;2 




(73) 



Finally, in [T^ the density profile in the point p = q {a, /?, 7, 6 arbitrary) was computed. This 
again corresponds to a simple representation of the algebra (see (pO|)). 

Much less is known about the partially asymmetric diffusion process. The current J in 
the large-L limit was determined in |jl5|. In analogy with the case of completely asymmetric 
diffusion discussed above, a phase diagram with only three different phases was proposed on 
the basis of the form of the current. The relevant variables for determining the phases are 
«;+(«, 7) and k+{P,S), where 
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K+{x,y) = —[-X + y + p- q + ^J{-x + y + p- qY + Axy] . (74) 
Note that for = ?/ = g, p = 1 this definition reduces to the one for K,j^{a) (see above). In 



terms of these variables the current phase-diagram exhibits the following three phases |15 
Phase A : k+(/5,5) > /?+(«, 7) , k^{P,6) > 1. In the limit L ^ 00 the current J is 
J=^^^^{{f3-S){p-q)-{[3 + 6)' + {f3 + 6)^{P-6-p + qy + Af36} . (75) 



Phase B : K+(a,7) > /t+(/3,(5) , K+(a,7) > 1. 
^ " 2[p- q) { ~ ~ ^) ~ + + (a + 7)\/("-7-P + g)^ + 4a7 } . (76) 



Phase C : k+(/3,5) < 1 ,/€+(«, 7) < 1. 

i = ^. (77) 

These results (which are the same as the corresponding mean-field results derived in Appendix 



D) are summarized in the phase-diagram shown in Fig. 2 |T5| 




1 2 3 4 5 



FIG. 2. 

Phase diagram for the current on a large lattice for p > q. The regions above the dotted 
line in phases A and B are accessible by finite-dimensional representations (see below). 
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We remark that in order to evaluate the current it is not necessary to first obtain a complete 
representation of (|54D, only certain matrix elements are required [0. As we are interested in 
general correlators we now turn to a detailed study of representations of (p^). 



VII. FINITE-DIMENSIONAL REPRESENTATIONS OF THE QUADRATIC 
ALGEBRA AND THE PHASE DIAGRAM 

As we have seen in the last section, the calculation of the density profile in the fully 
asymmetric case was done choosing the parameters of the problem such that one obtains an 
infinite-dimensional representation of a very simple form. We would like to do the calculation of 
the density and the correlation functions (which up to now are only known for a = (3 = p = 1, 
7 = 5 = g = 0, which corresponds to = K+(a,7) = in a large region of the 

parameter space. In order to do so, we will use the finite-dimensional representations of the 



algebra. What kind of results can one expect ? An inspection of the formulas (|59[) suggests that 
if one writes C = exp{—lH[), IM plays the role of a space-evolution operator although \ V) and 
{W\ are not eigenvectors of H. Two scenarios are possible. If C is diagonizable, one expects 
an exponential decay of the density and correlation functions. If C is not diagonizable one 
can anticipate an algebraic behaviour. From the mean-field analysis one expects an algebraic 
behaviour on the coexistence line and in the domain K_|_(a,7) < 1 and k^{/3,6) < 1. As we 
will see, the finite-dimensional representations will only access the coexistence line. Using 
finite-dimensional representations makes the calculation of the correlation functions relatively 
simple. However there is a price to pay, namely one has to solve the constraint equation which 
has a rather complicated expression in terms of the parameters a, j3, 7, 6,p, q. 

Let us first investigate the question which regions in the phase diagram given in Fig. 2 are 
accessible by finite dimensional representations. 

a) One-dimensional representation 

It exists whenever the constraint X7 = is fulfilled. In terms of the variables the 
constraint reads 

n^{(3,5) = ^^ (78) 

It is completely straightforward to evaluate the current and the density profile in this 
case. We find 

aP — 7(5 



a + p + -f + S 



^ a + p + -f + S l + /€+(«,7) l + K+{p,S)' ^ ' 

The second to last equality is established after some cumbersome computations using 
the constraint x-j = 0. 
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b) Two-dimensional representation 

Let us consider the 2-d representation in detail. The constraint /2 = is expressed in 
terms of the Xj's as 

{X2 - X3)[x7{xl - XiX^y + X2{X2X5 - XiXq){x2Xq - X4X5)] = . (80) 

Taking X2 ^ x^ {x2 = x^ corresponds to the unphysical situation p = —q) we can cast 
(PPP in the form of a quadratic equation for 7 

= C7^ + 67 + a , 

a = Pp'^(a'^(Pp + Sq + pq) + aq{Pp — 5p — p^ + 2Sq + pq) — 5(f{p — q)) 
b = pq{aj3p{l3 -5 -p + 2q) + a6q{(3 - 6 - 2p + q) - (36q{2p - q) 

—apq{p — g) — 5q^{5 + p — q)) 
c = -Sq'^(pp + qS + pq) . (81) 

which allows us to readily express 7 as a function of the other five parameters. Note 
however that the solutions of (pT)) still have to be supplemented by the condition 7 > 0, 
which excludes one of the two roots of (^) . The remaining one yields 

7i + 72 + 73 
7 = , 

74 



74 

71 = —p{—aP'^p + aP5p + ajSp^ — aPSq + aS'^q — 2aPpq + 2aSpq) , 

72 = —p{2l35pq + ap^q — a5q^ — [35q^ + 5'^q^ — apq^ + 5pq^ — 5q'^) , 

73 = p{al3p + a5q + apq + 6q^)sJ P'^ + 2^6 + 6'^ - 2(3p + 25p + p'^ + 2q{(3 - 5 - p) + q^ 
-fi = 2Sq{Pp + 6q + pq) . (82) 



For the special case 7 = 5 = (pOf ) has the simple solution P = —q + 

Using (|74D and (0) one can show that the region of the phase diagram accessible by the 
2-d representation is (without loss of generality we assume p> q) 

«+(/3,5)>— (83) 
K+(a,7) 

This can be easily checked for the case 7 = 5 = 0, for the general case we carried out a 
numerical analysis. The region described by (|83D covers the area above the dotted line 
in Fig. 2, i.e. most of the phases A and B. We note that for the case of symmetric diffu- 
sion p = q the two-dimensional representation does not exist. The infinite-dimensional 
representation is given by (^). 

c) Three-dimensional representation 

An analogous analysis of the constraint /s = for the three-dimensional representation 
leads to the same constraint (^). We believe this to hold for any finite-dimensional 
representation as well. 
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VIII. MATRIX ELEMENTS OF THE TWO-DIMENSIONAL REPRESENTATION 

Using (H) ; (PI) ; (|10|) cind (|14[) and performing a similarity transformation with 

S=[ ° I , (84) 
we obtain the following form for the matrices A and B 

A=\^^ ^ j , B=\^^ ^ j , l^) = I ^ 1 , (W^l = (l 0) , (85) 

where 

/i = — , 02 = — , a2 = — • (86) 

Using the constraint /2 = 0, 02 and 62 can be rewritten as e.g. 

{aP - j6){pP + q6) - (p - q)[pP{a + 6) + q5{[3 + 7)] 




a2 



p^aj3 — q'^'-)6 

{a/3 - 75) (pa + g7) - {p - q)[paif3 + 7) + i?7(a + 5)] 



(87) 



p'^aP — q'^'jS 
For actual computations properties of the matrix 

C = D + E= " + ^ + ^ + ^ + " + ^ A + -^±^B (88) 
a/5 — 75 a/5 — 75 a/5 — 7^ 

are of central importance. We have to distinguish between two cases: if a 7^ /5 or 7 7^ 5, C 
can be diagonalized 

SCS-' ={ ] , S = — 

^ _a + /5 + 7 + 5^ a + 7 ^ 

a/5 — -y6 a/5 — 7^ 
^ ^ a + /5 + 7 + (5 ^ /5 + ^ ^gg^ 

a/? — ■y6 a(3 — ■y6 

Note that the determinant det 5 = (/5 + 5)/i(A_ — A+) is different from zero unless A_ = A+ 
(the prefactor vanishes only if x-j = 0, in which case the 2-d representation breaks up into 1-d 
representations). In the case A+ = A_, C can no longer be diagonalized but only be brought 
to Jordan normal form. The condition A+ = A_ can be rewritten as 

(a + 7)02 - (/5 + 5)62 P-Q e / ,Q\ fx \^ fan\ 

= = pW-gV ~ " ~ • 

As the 2-d representation does not exist for p = g we conclude that the only case, in which 
A+ = A_ is if 

aS + p{a - (3) + q5 
' = W~q • ^''^ 



A numerical analysis of the two conditions ( ^21) and ( pi]) yields that a necessary condition for 
A+ = A- is that a = /5 and 7 = 5, i.e. the coexistence line of Fig.2. 
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IX. CORRELATION FUNCTIONS OFF THE COEXISTENCE LINE (A+ 7^ A ) 

The density profile (rj) is readily evaluated using (|59| ) ,(|88| ) and ( ^9]) 

{W\0-^DC^-^\V) {W\S-^SC^-^S-^SDS-^SC^-^S-^S\V) 



Af' \ ^ fXl'' 



{W\CL\V) {W\S-^SCLS-^S\V) 

(92) 



{W\ ^ \D\ \V) 





^ A:^ 

where \V) = ^QaW\ = ,,(,.,,)t(it,-^lf_(,^,),,i ( HP + 5)h , (« + 7)a, ), and 

n a+5+5b2 j 

Note that one of the matrix elements of D vanishes. This will have consequences for the shape 
of the density profile and the two-point function. 

The scalar products are easy to work out and lead to the following result for the density 
profile 

(tj) =Vt{uQ + ui exp {~~^^ + '^2 exp ^ ' (9^) 

where 

1 



[a + /5 + 7 + 5+(a + 7)02] [{13 + 5)1)2 - (a + 7)02 exp (^)] 
i = In i—\ = In ( '^ + P + ^ + ^ + iP + ^)bA 

^0 = i/3 + 5)b2{a + 6 + aa2) , 
uji = -a2&2 (a/5 - 7^) , 

uj2 = -{a + l)a2{a + S + 5b2) . (95) 
Two cases have to be distinguished: 

• C < 0, which corresponds to the case K+(a,7) > k+(/9,(5) i.e. Phase B. In this case the 
profile for L ^ is of the form 

(r,.) = m< + c< exp (^^^ + O (^exp (-^)) , (96) 

where 
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a + S + aa2 
a + i3 + 'y + 5 + {a + 7)02 



< [a + /5 + 7 + 5+(a + 7)a2](/3 + 5) ' ^ ^ 

The average density starts at the value m< at the left boundary, remains constant 
throughout the bulk, and eventually exhibits an exponential increase to the value m<+c< 
at the right boundary. 

C > 0, which corresponds to the case k+(q;,7) < 5) i.e. Phase A. In this case the 

profile for L ^ ^ is of the form 

0{exp(-^]] , (98) 




ICI 



where 



a + S + Sb2 
" a + (3 + -f + 6+{(3 + 6)b2 ' 

boiad — ■yd) , . 

c> = 7 o , ,^ 1, r < . 99 

[a + p + -f + 5 + {p + 5)b2]{a + -f) ^ ' 

Here the density starts at the value m> + c>exp(— ^) at the left boundary, increases 
exponentially to m>, and remains constant until the right boundary. 

Let us now take a closer look at the expressions for the bulk densities m< and m> and the 
correlation length \C\. It turns out that they are ^^universal" in the sense that they depend 
only on the two variables K+(a,7) and k,+ {(3,6) instead of all five independent parameters 
a, f3, 6,p, q. We start with m< and m>. One can show that 



m 



< 



1 + K+(a,7) 



1 + K+{(3,6) 

The equality of (|100|) with (|97D and (^) is established analytically only for = 7 = 5, and 
numerically to machine accuracy for the general case. As we will now argue, we believe (|100|) 
to hold not only for the 2-d representation but in general for phases A and B: in Appendix F 
a mean-field analysis of the partially asymmetric diffusion process is carried out, and m< and 
m> in phases B and A are determined. We denote the corresponding results (see Appendix 
D) by m^^MF and m^^MF- It is straightforward to demonstrate analytically that 

1 K^{j3,6) 

m<^MF = ^ 7 r , m>,MF = ^- rw^ ■ (101) 

l + K+(a,7) l + K+{p,d) 

This shows that the mean- field theory result is universal in phases A and B, and in addition 
is exact whenever the 2-d representation exists. Based on this observation and the results of 
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for = '-y = S = q (see eqn (0)) we conjecture that (|100|) holds true throughout phases 
A and B. 

The value of m> amd m< as a function of K+(a, 7) and S) is shown in Fig.3. Accord- 

ingly phase A is identified as a high-density phase with m> > |, and phase B as a low-density 
phase with m< < |. 



>1 o 





«:+(«, 7) 



FIG. 3. 

Plot of the bulk densities m>, m< in phases A and B in the region accessible by the 2-d 
representation . 

Unlike the quantities m> and m< the coefficients c< and c> of the exponentials are not 
universal in the sense that they are not only functions of n+{a, 7) and k+(/3, 6). However, the 
correlation length |C| in the exponential can be expressed as 

This shows that the correlation length \(\ diverges when «;+(«, 7) and n^{l3,6) approach the 
coexistence line: ( —>■ 00 when approaching the coexistence line from phase A, and ( — > —00 
from phase B. It is interesting to compare (^) with ( |102| ). Surprisingly enough, for phases Aj 
and Bj the = 7 = 5 = g correlation length has precisely the expression ( |102| ). This ceases to 
be the case for phases Ajj and Bjj. We would like to stress that although the mean- field values 
( |10CI| ) are exact, the correlation length ( given by ( p.02| ) cannot be obtained in the mean-field 



approximation. 

The two-point function can be evaluated in way analogous to the case of the one-point 
function discussed above. After some straightforward computations we obtain 
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{TjTk) = ^' (uj3 + uj4exp + ^sexp (- — ^ — ^^'j + c^eexp {~~^\] ' ^^^^^ 



where 



XlliP + 5)b2 - (a + 7)^2 exp (^)] 



^3 = (/3 + 6)1)2 



to'4 = —02^2 

= -a2&2 



a + 6 + aa2 
a[3 — ■y6 ' 

a + 6 + 51)2 
a[3 — 7(5 ' 



u^ = -{a + -i)a2[—^—^\ . (104) 

Again we have to distinguish between two cases 

• C < (low-density phase): The connected two-point function in the large L limit is given 
by 

{rjTk) - {Tj){Tk) = c<(m> - m<) exp (^—^^ " 4 exp [——^ — ' (105) 

It is different from zero only very close to the right boundary, from where it decays 
exponentially. 

• ( > (high-density phase) : 

{TjTk} - {Tj){Tk) = c>(m< - m>) exp I — - — I " c> exp I I , (106) 

Thus the connected two-point function in phase A is different from zero only very close 
to the left boundary, where it exhibits an exponential behaviour. 



X. CORRELATION FUNCTIONS ON THE COEXISTENCE LINE (A+ = A ) 

For the case a = /3, 7 = 5 we have k+(q;, 7) = fi;+(/5, 5) > 1 with 



/t+(a,7) = 




(107) 



and are thus on the phase boundary between the high-density phase A and the low-density 
phase B. Here we have taken into account the constraint (0) for the 2-d representation, which 
can be solved with the result 
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p 

J = -p + (a + q)J - . 

V Q 



(108) 



The matrix elements of the matrices A and B in 

02 = &2 = -ifl = 



P , 

1 — a 



which leads to the following form for the matrix C 

1 



C 



a — 7 



((2 + a2)/ + a2P) , P 



are found to be 

Vq + Vp 

-1 i 
i 1 



(109) 



(110) 



The matrix P has the property P^ = 0, which is important for carrying out the calculations 
below. Using the fact that = it is easy to show that 



(2 + a2 



'1 - k)a2 
ika2 



{a — 

which implies that the normalization is given by 



ika2 
2 + {l + k)a2, 



(111) 



{W\C^\V) = ^^^^^(2 + (1 - L)a2) 
[a — 7)-^ 

Using ( |112D we can easily evaluate the current 

J ^ / a-7 \ / 2 + (2-L)a2 \ 
V2 + a2/ V2 + (1-L)a2y ' 

In the thermodynamic limit this simplifies to 



(112) 



(113) 



J 



a — 7 

2 + 02 



pq- 



Q 



Q 



Q 



which in turn can be shown to be equal to (|76D . 
The density profile is 



{a + 7)(2 + ^2(2 - L)) + a^(7 - La) 
(a + 7)(2 + a2)(2 + (l-L)a2) 



+ 



P «;+(", 7) 



(a -7)^2 



(114) 



> + 7)(2 + a2)(2 + (l-L)a2)-' ' 



(115) 



In the thermodynamic limit j, L ^ (yo j- = x fixed, this turns into 



Q 



l + Q 



1-6 
l + Q 



X 



(116) 



This means that for p > q the density increases linearly from at the left end of the 
chain to y+q' whereas for g > p is decreases linearly from to y+q- Clearly the profile 
is symmetric under simulataneous interchange of p and q and left and right as it should be. 
The most remarkable feature of the profile ( |116| ) is the fact that it is independent of the 
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injection/extraction rate a. The only relevant parameter is the ratio | of the diffusion rates 
to the right and to the left. This fact is probably a feature of the particular representation 
we work with since in the completely asymmetric case the current and density profile on the 
coexistence line are dependent on the boundary condition a (see (0)), and we expect that 
in general the current and density profile in the partially asymmetric case will depend on the 
boundary conditions as well. 

The two-point function can be determined by using ( |111| ) in ( p9D and is found to be of the 
form 



where 

UJ-j 



(TjTk) =UJ7 + Usj + UJgk , (117) 

(a + 7)2[2 + 02(3 - L) + al] + 2a\Y{^ - La{a + 7)] + aH-f"^ - La^] 
(2 + a2)2(2 + a2(l-L))(a + 7)2 ' 
{a - 7)ai(a + 7 + 702) 



(2 + 02)2(2 + 02(1 -L))(a + 7)2 



(a-7)a2(a + 7 + aa2) 
(2 + a2)2(2 + a2(l-L))(« + 7)2 " ^ 

In the limit j, k, L ^ 00 with = x and j = y fixed this simplifies essentially 

{r..r,.> = ( —) + + -—^y . (119) 



Using (|116|) we finally arrive at the following result for the connected two-point function 

{TL.ny) - {TL.){TLy) = ([^) X{1 - y) . (120) 

Like the density profile the connected two-point function is independent of a. It takes its 
maximal value ( 2(1+^) ) ^^e middle of the chain, and decreases to zero at both boundaries. 
It is interesting to note that a similar expression for the connected two-point function was ob- 
tained for the completely asymmetric (particles hop only to the right) deterministic exclusion 
model with stochastic boundary effects 



XI. SUMMARY AND CONCLUSION 

In this paper we have first presented the Fock representations of the quadratic algebra. 
The representations can be either infinite-dimensional or finite-dimensional. Each finite- 
dimensional representation is characterized by a constraint on the seven parameters of the 
algebra. The matrix elements of the two generators of the algebra are given by recurrence re- 
lations. Only in the cases where these relations can be solved in a simple way the corresponding 
Fock-representation is useful for studying the physical problem of partially asymmetric dif- 
fusion with open boundaries. Quadratic algebras with more than two generators, which are 
relevant for many-state problems, will be presented elsewhere. The structure of the associative 
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algebra is different for these cases p8|. The boundary conditions define representations of a 
different type as compared to the Fock representations considered here. 

The quadratic algebra appears in the DEHP Ansatz for the steady-state probability dis- 
tributions of one- dimensional reaction diffusion problems with two-body rates and injection 
and extraction of particles at the ends of the chain. As shown in section V and also in P,P3| , 
the quadratic algebra appears also for some more general reaction-diffusion processes. As also 
shown in the present work, if one considers three-body rates, cubic algebras occur. We have 
not studied the representation theory for that case. 

We have used the 2-d representation of the quadratic algebra to compute the density profile 
and correlation functions. Both have a special dependence on the coordinates (see (|106D and 
(|105|) ). The parameter dependence is also peculiar. Certain quantities like the density around 
the middle of the chain or the correlation length depend on the parameters a, j3,'y,6,p,q 
only through two functions /€+(a,7) and 5) defined in the text. Our results and those 

of [|ll|,|l^,|l^] suggest that the density in the bulk is given by the mean-field prediction in the 
thermodynamic limit. Again from our results and those of [^,^,0 it looks like the correlation 
lengths again are "universal" (they depend only on K^{a,'-f) and K,_^.{j3,6)) in regions Aj and 
Bj of the phase diagram. 

By means of (^91) our results ( |105| ), (|106|) and ( |12CI|) for the connected two-point functions 
can be directly translated to the XXZ quantum-spin chain (|47|) with non-diagonal boundary 
terms. This yields new and nontrivial results for correlators in the XXZ chain with boundaries. 
We have not considered the problem of time-dependent correlation functions. Some recent 
numerical results of Bilstein |2^ can give a hint in this direction. He studied the finite-size 
scaling bahaviour of the energy gaps of the quantum hamiltonian. It turns out that in the low- 
and high-density phases the system is massive. In phase C (see Fig. 2) the system is massless. 

3 

Both real and imaginary parts of the energy gaps vanish like L~2, where L is the size of the 
system. The coexistence line is also massless but features a less simple length dependence. 

Finally, as a by-product of our work on the DEHP construction, we have given a simple 
way to construct irreducible representations of the Uq{SU{2)) quantum group. 
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Appendix A 

We would like to show that in the Fock representation of the quadratic algebra (0) the 
matrices A and B can be written in tridiagonal form. This is a property of the Fock represen- 
tation only, and is independent of the parameters Xj. Instead of presenting the general proof, 
we give two examples which will help the reader to see the mechanism behind the proof. We 
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start with the 2-d representation. We choose \V) 

Wi such that (W^|V^) 7^ 0. We then perform a similarity transformation 
1 a 



and take {W\ = {wi, W2) with 



S 



1 



^ S\V) = \V) , {W\S-^ = {wr, -awr + W2) . 



(121) 



where we choose a = ^ and Wi = 1, which results in {W\ = ( 1, 0) and = 1. In this 

basis A and B now must be of the form 



A 



B 



ai2 

We now perform another similarity transformation 

1 





^21 ^22 . 



(122) 



S 



.0 (3, 



{W\S~^ = {W\ , S\V) = \V) 



(123) 



where we choose 



1^ and find two equivalent representations 



A 



±Vai2&2i 
a22 











B 



, i V '3'12&21 ^22, 



(124) 



Here we have assumed that both au and 621 are different from zero. One can show that if 
ai2 = or 621 = the algebra (H) yields the condition = 0, which corresponds to the 1-d 
representation and thus the 2-d representation is not interesting. Renaming the entries of A 
and B we arrive at 



A 



/i 
02, 



B 







(125) 



Let us now consider the 3-d representation. We can repeat the similarity transformations 
from the 2-d case to bring A and B to the form {f[ 7^ 0) 





(0 f[ 


9[\ 




(0 





\ 


A = 


a'22 


«23 


, B = 


fi 


b'22 


b'23 




[0 a',2 






\9[ 


b'32 





(126) 



where \V) 



/1\ 




and {W\ 



0, 0). Taking the further similarity transformation 



/I 



S 







\ 



cos(e) sin(0) 
Vo -sin(^) cos(^) y 
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, tan(6') 



9[ 
fi 



(127) 



we get 



A 



/O /i \ 

0,22 '3'23 
V 032 033 / 



B 



/ \ 

/l ^22 ^23 
V 632 &33 / 



128) 



Finally, by means of a diagonal similarity transformation we can bring A and B to the form 

/0/iO\ /000\ 



A 



a2 /2 
VO h2 03/ 



B 



fi b2 k2 

V /2 ^^3 7 



(129) 



This procedure can be generalized to any dimension of the representation. We can thus 
search for representations of the quadratic algebra starting with the tridiagonal form 

/ 


A= h2 as f3 0... (130) 


V... 



/l 














o...\ 


a2 


/2 











0... 


h2 


03 


fs 








0... 







04 


h 





0... 



/O oo...\ 



B 



fi b2 

/2 




V... 



ko 0. 



/3 



0. 
k^ 0. 



(131) 



Appendix B 

In this appendix we discuss a generalization of the DEHP-formalism to chemical processes 
of the types given in (|35D that incorporate three neighbouring sites. By construction all 
models discussed above are contained in the present formulation. From a physical point of 
view this generalization is quite interesting. For purely diffusive processes it allows us e.g. to 
let particles hop to the next-nearest neighbour site if the nearest neighbour site is unoccupied. 
In the corresponding traffic-flow picture this corresponds to letting cars move faster or slower 
depending on whether the road ahead is free for a long or a short distance. 

We consider a master-equation of the form 
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dP ^"^ 

= ^ = - 51 H {Hk^k+l,k+2)ll'Jlll'Jlll Pdn, T2... Tfc^i, 7fc, 7fc+i, 7fc+2, rfc+3 . . . Ti) 
fc=l 7fc.7fc+l.7fc+2 

-E(M;7.^^L(7i72r3...r,.)- E (/^l-il)?^-,^^ PlI^i . . . r^^sTL-m) , (132) 

7i72 7l-i7l 

where 

( rr ^7.,7.+l,7.+2 _ / ^^'1+1 '1+2 7, = T", , J = fc, A; + 1, A; + 2 

(7ifc,fc+i,fc+2j^ - < {166} 

k, k+2 I _r7fc.7fc+i.7fc+2 p op 

Here r°f^,^, is the probabihty per unit time that the configuration (a/37) three neighbouring 
sites changes to the configuration (a'/3'7') . Following the analysis in section 1 above we demand 
for the probability distribution in the bulk that 

E (^k,k+l,k+2)ll' Jill' Jill Pl{Ti, T2... Tk-l, -fk, 7fe+l, 7fc+2, Tk+S ■ ■ ■ Tl) 
7fei7*;+l)7fe+2 

= Xr^Pi-liri . . . Tk-lTk+l ■■ - Tl) - Xr^^^Pi-lin . . . TkTk+2 ■ ■ ■ Tl) 

+yr^PL-l{Ti . . . Tk-lTk+1 ...Tl) - yr^^^PL-l{Ti . . . Tfc+iTfc+s . . . Tl) 
+Zr^^^PL_l{Ti . . . TkTk+2 ■ ■ ■ Tl) - ^^..^^Pi-i (Ti . . . Tk+lU+s • • • ^l) 

+tr,r^+,PL-2{ri . . . Tfc^iTfc+a ■ ■ ■ Tl) - tr,^,r,+^PL-2{Ti . . . TkTk+S ■ ■ ■ Tl) . (134) 

This ensures that there will be no contribution from the bulk to the r.h.s. of ( |132| ), as all 
terms will cancel when summed over {k,k + l,k + 2). As compared to (|53D it is now possible 
to include terms containing probability distributions Pl-2 with two fewer particles. The r.h.s. 
of ( p.34| ) can be simplified by introducing the notation 

/^bfc = + Vbk ' = Zb^ - Xb^: , Afc^ = -/ifo^ - Z/fe^^ . (135) 

From the boundaries we get the conditions that 

E (^12)nr2 ^i(7l72T-3 ■■■Tl) = U^r2PL~2{n ■ ■ ■ Tl) + /irii^L-l (^2 ■ ■ ■ Tl) 
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-\r2PL-l{riT'i ...Tl) 
E {hL-lL)1lzJltPL{r, . . . TL-21L-11L) = 

7L-17-L 

-tr^_irLPL-2{Ti . . . TL-2) " /iri,_i -Pl-1 (^l . . . Tl-2Tl) + KlPl-i{ti . . . Tl-i) • (136) 

Here we have allowed for arbitrary processes of the type introduced in (^) to occur on the first 
two and last two sites of the lattice. Putting everything together we arrive at the following 
algebra 
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-n 



( ) 




( 





\ 


D^E 




fiiDE + viDE + Ao^*^ + tiiE 


-hoD 


DED 




fiiED + uoD^ H 


VXiDE + ti^D 


-toiD 


ED^ 




fioD"^ + viED - 


h XiED + tQiD 


-tuE 


DE^ 




fiiE^ + uqDE H 


h \qDE + ti^E 


-tooD 


EDE 




^qDE + viE'^ - 


V XqED + toiE 


— tioE 


E^D 




^loED + uoED + AiE2 + tooD 


— toiE 






\ 





J 
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where 



n 
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pOOl 
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pOOl 
000 



pOOO \ 

J- 111 \ 
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pOOO 

oil 
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rrOOO / 
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:i38) 



It is clear that only seven of the eight equations in ( |138| ) are independent. The boundary- 
conditions are rewritten as 



12. 



TlT2 



{W\[^iD + (1 - 7i)E] + (1 - 72)E] 
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{W\ 



firAr2D + (1 - r2)E] - K,[riD + (1 - n)E] 



^ ihL-lL)l',-_rr^bL-lD + (1 - ^l^,)E][^lD + (1 - ^l)E]\V) 



- t 



TL-lTL 



+ XrA^L-lD + (1 - rL-l)E] - fir^^A^LD + (1 - Tl)E] 



\V) 



(139) 



Note that the algebra ( [L38| ) is cubic and that the conditions at the boundaries are quadratic 
in D and E. 

If we constrain ourselves to diffusion processes only, the system ( |138| ) decouples into two 
sets of respectively two independent equations 

/ D^E \ 

-XiDE + XqD^ + tuE - tioD 



Hi 



DED 
\ ED^ J 



fXiED + uqD^ + XiDE + ti^D - t^iD 



31 



where 



/ DE^ \ 

EDE 
\ E^D ) 



-Hi 



/iiE^ - ^^DE + UqE - tooD 
^HoDE + uiE"^ + XqED + toiE - ti^E , 
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pllO _ pllO 

101 oil 



piio 
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pioo 

010 



pioo 

001 



pioo 

010 



rioi 

110 

_pioi _ pioi 
^ 110 oil 

pOlO 
100 

_p010 _ poio 
100 001 



pOll 
110 

pOll 

101 , 
pOOl ■ 
100 

pOOl 
010 . 



(141) 



In order to demonstrate that there exist solutions to these equations we will prove the 
existence of a one- dimensional representation for the special choice of boundary conditions 



—a 



(142) 



which correspond to injection of particles at sites 1 and 2 with probability a if both sites are 
empty, and extraction of particles at sites L — 1 and L with probability (3 if both sites are 
occupied. The one-dimensional representation exists ii D = E = 1 and 



a 

pllO I pOll 
101 "T 101 



pllO 
101 



P 

plOl I plOl 
oil "1" 110 ' 



pllO 

oil 



plOl _ pi' 

-■-110 1 
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plOO I pOOl _ pOlO I pOlO 
010 "T 010 — nm "r J- 
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010 ,n 
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pOlO 
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The corresponding probability distribution is trivial 



(143) 



(144) 



which means that all configurations are equally represented. Although the existence of the 1-d 
representations is a nontrivial fact, the corresponding physics is not particularly interesting. 
It would be very interesting to construct finite or infinite dimensional representation and use 
them to compute currents, density profiles etc. 

Appendix C 



In this appendix we consider the special case = 7 = 5 and a = (3 = p — qoi the algebra 
([5^). As was shown in (^) the quadratic algebra reduces to a Q-osciallator algebra. We 
introduce the following notation for Q-numbers 



and Q-binomials 



Cl{Q) 



{n] 



M{n-1}...{2} 



l-Q 



n: 



" i(i 



{n - p]{n -p-l]... {2}{p]{p - 1} ... {2} 



m=i(i 



Q^)nr=f(i-Q* 



(145) 



(146) 
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In order to compute the current or correlation functions a convenient representation of 
{D + E)^ is required. In terms of a and we find 



2 _^ a + aT 



N 



N N 



\1-Q ^/Y^, 



^C^(l)(a + a 



Wo 



(147) 



In order to evaluate e.g. the normalization (0|C^|0) it is convenient to decompose powers of 
a + into normal-ordered expressions defined via 



p=0 



The decomposition is of the form 



:i48) 



[a + = E ^n"^ ■ (« + 



^\n—2m . 



m=0 



Using the identity 

(a + a^) : (a + a^)" :=: (a + a^)"+^ : +M : (a + a^)""^ : , 
one readily obtains recursive expressions for M^'^^ 

n— 1 n— 3 n— 5 

1=1 1=1 1=1 

We first note the result for the cases Q = 1 (no deformation) and (5 = 

Mi'") 

n — 1 



n \ rn-l 

2m / fc=o 



(n-2m)! (2m)!! 

72 — 1 

Q=o V m j \m — 2 
After some tedious computations we find for arbitrary Q 



Mi™) 



1-g 



Mi^) 



1 



^ I "LA _ _ 2}) - {n} ( i^V^ - - 3} 



{2} 



(149) 



(150) 



(151) 



(152) 



(153) 



We did not succeed in obtaining a closed form for M^^^ for general k. The difficulty can be 
traced back to the recurrence relation ( |151| ) which although it is written entirely in terms of 
Q-numbers does not have Q- number solutions (the sum of Q- numbers is not a Q- number). 

Appendix D 

In this appendix we would like to show that the DEHP-Ansatz gives a new way to construct 
irreducible representations of the quantum group Uq{SU{2)) and that the quantum plane 
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appears in a natural way. We consider the special case = ^ = 77, A = g/j9, Z2 = vKv ~ <?) 
of (^41). The boundary conditions are taken such that = 7 = 5, then the limit a ^ 0, 
/? — >^ is performed. For = q; = /? = 7 = (5 the XXZ quantum spin hamiltonian ( ^T] ) 
corresponding to the diffusion process is invariant under the quantum algebra Uq{SU{2)) ( 
[P!7|,|5| and references therein). The stationary state is no more unique since the ground 
state of the ferromagnetic chain is L + 1-times degenerate corresponding to a multiplet of 
slq(2). As the boundary conditions (Q) become ill-defined in the case = a = /3 = 7 = (5we 
carefully take the limit a — >^ 0, /5 — > with = 7 = 5 as follows 

a = -, (3 = % D = eD, E = eE, 
e a 

{W\E = e{W\, D\V) = d\V), (154) 
where D and E are seen to obey the quantum plane equation 

DE = Q^ED , Q2 = - . (155) 

P 

The quantum-mechanical state corresponding to the stationary state of the diffusion process 
is obtained by applying (|38|) 

|0) = ^Erf^-V E Q2(-+-+-+^-^)nar| TT ... T), (156) 

fc=0 h<i2<...<ik 1=1 

where is a normalization factor. Note that |0) is a linear combination of L + 1 independent 
wave functions. In order to get the ground state of the XXZ chain (^) we still have to perform 
the similarity transformation (^), which changes the probability distribution of the stationary 
state (BII) into 



{W\ n {r,AQ'-'D + (1 - t,)e) \V) = A^Q'^iWl f[ (^r^D + ^^^E^ \V) . (157) 

This yields the following result for the similarity-transformed state \0)u, which is the ground 
state of the f/q(S't/(2))-invariant XXZ hamiltonian { ^7\ ) 

1 L k 

\0)u = —A'-Q^ E A^'t^'^^'e'^ E 0!^-^'^+-+'^-^' n ^.^ I TT • • • T)- (158) 

■'^ k=Q ii<i2<-<ik l=i 

We note that the k'th term in the sum is proportional to the state (>S'~)'^| TT • • •) obtained by 
acting with the quantum-group generators 



1=1 



We believe that using the DEHP Ansatz in order to get irreducible representations of quantum 
groups may have other applications. 
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Appendix E 



Here we address the problem of existence of solutions to the system of relations 
(|33|),(^),(|6^)-(|6^). It is convenient to define symmetric and antisymmetric combinations 
of rates 

p± _ pOl I plO p± _ pOl _|_ plO p± _ pOl _|_ plO 
a ~ 11 =■= 11 ) 6 ~ -"^ 00 00 ' ^ c ~ J- 10 01 ' 

■pi -pii I -pii -pi pOO _|_ pOO A Oi + ^ /'1«n^ 
= Tio ± , r, = ± , A = —— . (160) 

P + 7 

Note that due to positivity of the rates all symmetric combinations are automatically positive. 
With 

aP — 7(5 ' aP — ■y6 ^ ^ 

it is straightforward to show that 1-d representations exist under the condition that 

r+ + r+ = r+A + r+A-i, 
A(rJ^ + r+) = r?;A-i + C 

r;A-^ + Ar^ + 2r; + r- + r,- = + a-^) . (i62) 



It is obvious that one can choose the 12 rates such that equations (|162|) are satisfied. 

Appendix F: Mean-Field Analysis 

In this appendix we give a summary of the mean-field analysis for the partially asymmetric 
diffusion process on a lattice with L sites. Our discussion follows [|ll|)i3l' which deals with 
the completely asymmetric case. The particles hop with rate p (q) to the right (left) and are 
injected (extracted) with rates a (7) at site 1 and rate 6 (/3) at site L. 

In a stationary state the density at site j is time-independent, which implies 

= ^ = (g - P){rj^irj) + {p~ q){rn+i) + P{r,-i) + g(r,+i) - [p + q){T,) . (163) 

Denoting (rj) by tj and decoupling the two-point functions (rjTfc) = tjtk leads to the following 
set of mean-field equations 

ptj-i + qtj+i - ptj^itj - qtjtj+i = {p + q)tj - ptjtj+i - qtjtj^i , j = 2 . . . L - 1 (164) 
a(l -ti) + gt2(l -ti) =7^1+^^1(1-^2) , 

ptL + gtL(l - tL-l) = S{1 - ti) + ptL-i{l - tL) . (165) 



The bulk equations ( |164| ) can be rewritten as 



tj+itj = —tj+i H ^^tj + c , (166) 

p - q p - q 

where c is an arbitrary constant (related to the current J). The net mean- field current from 
site j to site j + 1 is defined as 
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J = ptj{l - tj+i) - q{l - tj)tj+i = -{p - q)c , (167) 
and is independent of position as it should be for a stationary state. If we define 

,. = I^L + JL\ ^ (168) 

^ p + q\^ p-q) ^ ' 

the new quantities Sj are seen to obey the recursion 

d , pq (p-qV 

Sj [p + qr \p + qj 

_ pa + q^y + pq - {p + q^c' _ p6 + q(3 - pq + {p + q)'^c' 

" (p + g)(« + 7) ' " iP + q)iP + S) • ^ 

We note that < c' as otherwise S2 > 1, which is unphysical as < tj < 1 , Vj (which imphes 
^ Sj < ^H^)- -'^^^^ completely asymmetric case we now have to distinguish three 

cases: 

. c'>\ 

In this case the recursion has no fixed point and Sj would eventually turn negative, which 
leads to an unphysical solution. Thus we can exclude this case. 

This case corresponds to the maximal current phase C in the phase diagram (see Fig. 2). 
The current is given by (see above) 

J=-(p-,)c=^.'-^ = ^, (170) 

[p — q) P — Q 4 

which is the same as the exact result ([77|). 



In this case the recursion (|169|) has two fixed points 

1 



s± = -[l±VT^) ■ (171) 

Writing Sj = Vd^^^ we see that the u/s are subject to the recursion 

1 , , 

Uj+i + Uj^i = -j=Uj , (172) 

which is recognized as a special case of the recursion relation for Chebyshev polynomials. 
Eqn ( |172D is solved formally as 

M„(^) = sin[(n- 1)0 + 0] ,^ = arccos (^^) , (173) 
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(note that 9 is complex) which leads to the following expression for s„, j = n . . . L 

sm[n6 + (f)] 



sin[(n - 1)6* + 



(174) 



Using the two boundary conditions ( |169|) in ( |174|) completely fixes the values of and c' 
as functions of a, P, 7, S, p and q. For simplicity we introduce the notation si = di — d2c' 
and sl = + d^c' . The boundary condition for si implies that 

2{2d,-l)cos\e)-d2 
'''^'^^ = ^) ' ^^^^^ 

whereas the one for sl yields 

, . ^ 2 cos(g) cos(Lg) - (4 + 44 cos^(g)) cos[(L - 1)9] 

^'^^ 2cos(e)sin(L^)-(d4 + 4d3Cos2(e))sin[(L-l)e] ■ ^ ' 



Equating ( |175| ) and ( |176| ) we obtain 



= di(l - 4) sin[(L + 3)^] + [{2di - c/2 - 1)(1 - 4) + 2c/i(l - 2^ - 6/4)] sin[(L + 1)^] 
+ [{2di -d2- 1)(1 - 2^3 - ^4) + (rfi - 1)(1 - rfs) - dids] sin[(L - 1)^] 
-[(2rfi -d2- 1)4 + (1 - rfi)(l - 24 - (^4)] sin[(L - 3)9] 

+d3{l - di) sin[(L - 5)^] . (177) 

In the large-L limit (|177D turns into a fourth order polynomial equation in 2; = exp(2i^) 
{9 is complex). The polynomial equation can then be solved explicitly for 9 (and thus 
c') as a function of a, P, 7, 6, p, q. 



However, there exists a much simpler way to determine the mean field current |jTT[: in 
the low-density phase B we start out infinitesimally close to the unstable fixpoint s_, i.e. 
Sl = S- + e. The density stays at s_ throughout the bulk and only deviates towards s+ 
at the right end of the chain. Using the fact that c' = s_(l — s_) in the expression (|169|) 
for Sl, and then setting si = s_ immediately yields s_, and thus also c', as a function of 
p, q, a and 7. Inserting the resulting expression for c' into ( |170| ) then yields the current 
as a function of a, 7, p and q. The result found to be identical to the exact expression 

m- 

An analogous analysis can be carried out in the high-density phase A. Again the result 
is identical to the exact expression ( [75| ) 

We also can use mean- field theory to determine the density profile in phases A and B. From 
our discussion above it is clear that in phase B the density profile in the bulk is essentially 
constant and equal to the value of the unstable fixed point (we switch back from the Sj-variables 
to tj variables) 

«. = i(l-yrT4^) = i(l-jl-i^). (178) 
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Analogously, in phase A the profile in the bulk is constant and equal to the value at the stable 
fixed point t+ 



U = ]-{1 + Jl-^). (179) 
2 V P ~ Q 

Using the expression for the currents these values can be determined explicitly 



1 1 (p-qY 

t- =: m<,MF = \ a{a + j)K+{a, 7) + 7[a + 7 + p - g] , 

2 p — q \ 4 

t+ =: m>,MF = 1 + _ ^(^ + s)k+{P, 6) + 6[(3 + 6 + p-q] . (180) 

These expressions coincide with ( |100| ). Finally we note that the mean-field result for the 
correlation length ( does not repreduce ( |102|) and is incorrect. 
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